Dark Matter, Modified Gravity and the Mass of the Neutrino. 
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It has been suggested that Einstein's theory of General Relativity can be modified to accomodate mismatches 
between the gravitational field and luminous matter on a wide range of scales. Covariant theories of modified 
gravity generically predict the existence of extra degrees of freedom which may be interpreted as dark matter. 
We study a subclass of these theories where the overall energy density in these extra degrees of freedom is 
subdominant relative to the baryon density and show that they favour the presence of massive neutrinos. For 
some specific cases (such as a flat Universes with a cosmological constant) one finds a conservative lower bound 
on the neutrinos mass of m v > 0.31 eV. 



I. INTRODUCTION 

There is compelling evidence that the baryons in the Universe 
are unable to generate the gravitational potentials that we ob- 
serve on a wide range of scales. A simple paradigm can be 
used to explain this mismatch between light and gravity: the 
Universe is filled with an appreciable amount of matter which 
is cold (i.e. has non-relativistic velocities today) and does not 
interact with light. It has been shown that Cold Dark Matter 
(CDM) can explain a host of observation, from dynamics of 
clusters to the formation of the cosmic web JH . 

The CDM paradigm has been proposed within the context 
of Newtonian gravity and Einstein's theory of General Rela- 
tivity. It has been argued that these theories may not be valid 
on all scales. Indeed, proposals for modifying gravity have 
been shown to fit much of the currently available data 0]. 
A plethora of covariant theories have been studied in detail; 
TeVeS gravity, modified Einstein- Aether theories, conformal 
gravity, higher derivative actions, etc, have been advocated as 
possible rival theories to the CDM scenario |S|j,|J,|aLZll8l|9t]- 
There has been considerable effort in studying the cosmologi- 
cal consequences of these theories ITol [Til [l2l [l3ll . Given the 
level of precision of current cosmological data, it is possible to 
find severe constraints on these alternative theories and com- 
pare their ability to describe nature with the CDM scenario. 

There is an important, generic feature of covariant theo- 
ries of modified gravity which is often overlooked: although 
they tamper with the gravitational sector of the equations of 
motion, they also inevitably lead to the introduction of ex- 
tra degrees of freedom which may be interpreted as an exotic 
form of dark matter. Let us exemplify. Theories which mod- 
ify the Einstein-Hilbert action by, for example, replacing the 
Ricci scalar, R, by a function of different curvature invari- 
ants, f(R, R a pR al3 , • • • ), introduce higher derivative terms, 
and hence new modes. These new modes will contribute to the 
overall energy density. This is patently obvious in the case of 
theories where / is simply a function of R; such theories can 
be mapped onto normal Einstein gravity with an additional 
scalar field. This also true of conformal gravity, where the ac- 
tion is now constructed from the Weyl tensor. A field must 
be added to fix the scale of gravity and the resulting low en- 
ergy equations are fourth order [14]. More modern attempts at 
constructing theories of modified gravity have the same char- 
acteristics in a much more explicit way. In TeVeS [5], a scalar 



field and a vector field is introduced which not only mod- 
ify the gravitational field equations but also source the very 
same field through their stress energy tensor. In generalized 
Einstein- Aether theories, a time-like vector field is introduced 

a. 

Given what we have just said, there is an obvious ques- 
tion: aren't these extra degrees of freedom simply a contrived 
form of dark matter? It is conceivable that the extra degrees 
of freedom in modified theories of gravity may play such a 
role. If so, dark matter has been introduced through the back 
door. It turns out that the role of extra degrees of freedom 
in theories of modified gravity is more complicated than one 
might expect. In Skordis et al HlOfl . it was shown that the extra 
degrees of freedom in TeVeS can make a negligible contribu- 
tion to the background (or overall) energy density. Indeed, if 
TeVeS is to be consistent with big bang nucleosynthesis, the 
fractional energy density in these extra degrees of freedom, 
fix, must be under a percent. Yet even though fix <C 1, fluc- 
tuations in the extra fields could have a significant impact on 
the growth of structure. In particular, due to the modified na- 
ture of gravity, they could source the growth of gravitational 
potentials and sustain them through Silk damping at recom- 
bination. These results were corroborated in Dodelson and 
Liguori 11511 . where the fluctuations in the vector field were 
found to play an important role. 

Hence some theories of modified gravity can fit current ob- 
servations of large scale structure, either from galaxy surveys 
or the cosmic microwave background, even though fix <C 1. 
We would like to point out that the latter property is not 
generic. In some incarnations flB -C fix — 1 where fls 
is the fractional energy density in baryons. These theories end 
up being a hybrid of the two paradigms, modified gravity and 
dark matter, and in principle should be harder to distinguish 
from dark matter theories (although there are some sugges- 
tions of specific tests) QlMl- 

In this paper we will try to expand on an important fea- 
ture of TeVeS pointed out in Skordis et al : if one assumes 
that the Universe is flat and the only form of non relativis- 
tic matter consists of 5% baryons (consistent with Big Bang 
Nucleosynthesis), the angular power spectrum of the Cosmic 
Microwave Background (CMB) will differ significantly from 
observations. The only way to resolve this discrepancy is to 
introduce some form of non-relativistic matter, and the only 
one allowed within the known menagerie of fundamental con- 
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stituents of the Universe is a massive neutrino. To match 
observations of the CMB, neutrinos with a mass of approxi- 
mately 2 eV are needed. This result is clearly a hint of a more 
general statement that may be made about theories of modified 
gravity in which the extra degrees of freedom play a subdomi- 
nant role: if these theoriess are to agree with measurements of 
the CMB then they require the presence of massive neutrinos. 
We wish to see if this implies a lower bound on the mass of 
the neutrino. 



H. AN APPROXIMATE THEORY AND COSMOLOGICAL 
OBSERVABLES 

Let us consider a generic modified theory of gravity in the 
limit of homogeneity and isotropy. The physical metric (i.e. 
the metric which is minimally coupled to the matter fields) 
can be parametrized in terms of a scale factor, a(t) which has 
a logarithmic derivative, H = d In (a) / dt. The energy density 
of the Universe can be split into the normal degrees of free- 
dom, p (such as baryons, photons, neutrinos and dark energy) 
and the extra degrees of freedom that arise from the modifica- 
tions, px ■ The modified Friedman equations look somewhat 
like 



F(a,H)H 2 



8ttG 



(p + Px) 



(1) 



where G is Newton's constant and FH 2 is a function that 
arises from varying the action for a particular theory. In fact it 
is convenient to rewrite the equation in a more familiar form 
by defining an effective Newton's constant G e // = G/F. For 
the purpose of what follows we use a parametrization such 
that G e ff ~ G e ff(a); with a sufficiently flexible choice of 
parameters we can encompass cases where G e f / depends on 
a, H, etc. 

We consider a sub class of the theories, in which px < Pb, 
(where pg is the Baryon density). We consider a parametriza- 
tion such that 



PX fBPB + fRPR 



(2) 



where p^ is the energy density in radiation. We have that 
fs < 1 and the correct abundance of light elements requires 
that fn < 10 -2 . We also include in p, a component that be- 
haves like dark energy, pDE, with an equation of state w < 0. 
We find it convenient to parametrize the equation of state of 
the dark energy component as w = wq + W\z/(1 + z). Modi- 
fications to the gravitational sector may lead to accelerated ex- 
pansion at late times (such as those proposed in ffll^l8l [l2ll ). 
meaning that the dark energy could also arise from the extra 
fields in the modified gravity sector. Our dark energy term 
includes all of these possibilities. 

With the evolution of the scale factor in hand, there are 
a few observables that we may now calculate. Let us start 
off with the position of the first peak of the angular power 
spectrum of the Cosmic Microwave Background (CMB). It is 
a direct measure of the angular diameter distance and hence 
of the expansion history of the Universe from recombination 



until today and is the centrepiece of the analysis of this pa- 
per. Schematically we have the following picture II 1611 . Be- 
fore recombination (which occured at time t*), photons and 
baryons were tightly coupled and underwent acoustic oscilla- 
tions. During tight coupling the photon density contrast in the 
conformal Newtonian gauge, obeys the differential equation 
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where 



2>pb + 4/9 7 a 



4p 7 



(3) 



—— — —j is the sound speed, S[&, $, is a 
source (a function of the gravitational potentials, $ and ty) 
and derivatives are with conformal time r. In the WKB ap- 
proximation 01611 the two linearly independent solutions to the 
homogeneous part are cos kr s and sin kr s , and depend on the 
sound horizon r s (r) = f Q c s dr. The important thing is that 
these solutions are valid for any theory of gravity for which 
photons and baryons see the same physical metric. All modi- 
fications are to gravity or additional fields and implicitely alter 
the inhomogenous part through the source term S which only 
depends on the potentials $ and of that same physical met- 
ric. 

The physical scale, d* of the acoustic waves is set by 
the sound horizon at time t*, i.e. d* ~ J Q * c s (t)dt. Af- 
ter recombination, photons decoupled from the baryons and 
freestreamed towards us, travelling a distance given by do — 
ca o It° dt/a — da{a 2 H) where the subscript labels to- 
day. The angular size on the sky of the sound horizon, 6*, is 
given by 9* ~ !!a! g L The sound horizon at last scattering leaves 
a very distinct signature on the angular power spectrum of the 
CMB: a series of peaks and troughs. The spectrum generated 
at recombination is related to the spectrum today via a projec- 
tion through a spherical Bessel function je(k(r — r*) , where 
adr = dt, (an ultraspherical bessel function in the curved 
case). Once again this is independent of the theory of gravity. 
The position of the peaks and troughs in the angular power- 
spectrum is primarily dependent on the cosmological shift pa- 
rameter, 1Z which is related to the angular diameter distance 
and is given by 



n = - 



1 tt M (a ) 



(ao) 
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y = s/\^k{oo)\ 



da 



J (Ei«i(a)) 1/2 
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where f2, (a) is the fractional energy density of component i 
as a function of scale factor (K corresponds to the curvature 
and M to non relativistc matter) lfl7ll . 

Very few assumptions have gone into this calculation: the 
Universe underwent recombination and the horizon structure 
is a result of the expansion rate of the Universe. Current mea- 
surements of the CMB have reached a level of precision such 
that it is practically impossible to deviate from this simple pic- 
ture II X 911 . Attempts at changing these fundamental assump- 
tions inevitably lead to radical departures from this simple pic- 
ture and a gross mismatch to the data. So any theory of mod- 
ified gravity must lead to the basic picture of the CMB that 
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Model 


68% CL 


95% CL 


ACDM 


0.58 < m v < 1.17 


0.31 < ra v < 1.48 


wCDM 


1.58 < m v < 2.59 


0.92 < ra v < 3.02 


wCDM + Q K 


0.07 < m v < 1.05 


m v < 1.8 


ACDM+G(z) 


0.018 < m v < 0.62 


m v < 1.04 


wCDM + G{z) 


1.05 < m v < 2.33 


0.22 < ra v < 2.68 


wCDM +Q K + G(z) 


ra„ < 0.62 


ra v < 1.32 


ACDM +G(a,7,z) 


0.034 < m v < 0.71 


m„ < 1.05 



TABLE I: Results for different cosmological models for a compila- 
tion data set. 



we infer from the data. Hence we can use our Eq. Q]to work 
out 1Z for theories of modified gravity with the caveat that the 
fractional energy densities must be rescaled by the effective 
Newton's constant, i.e. we must replace fli by (G e ff /Go)ili 
in Equation [4] Throughout this analysis we consider a con- 
servative bound on the shift parameter: 1.63 < 7Z < 1.76 

Another useful observable, as measured from the Hubble 
diagram of distant supernovae, is the luminosity distance, d^- 
It is related to the angular diameter distance, (La, described 
above, through (Ll = (1 + z)dA, where 1 + z — ao/a defines 
the redshift at a given value of the scale factor. While the 
CMB gives us one measure of cLa at z ~ 1100, the Hubble di- 
agram of distant supernovae gives us a series of measurements 
of d,L out to z ~ 1.8. Again, as above we can use our modified 
Friedman equations to calculate and compare to the data. 
Given that we do not have to use any information about pertur- 
bations about the background, we make even fewer assump- 
tions. We use the group of supernovae, termed the 'gold' set, 
from the HST/GOODS programme BOh . complemented by 
the recently discovered higher redshift supernovae, reported 
in iH. 

Finally, we consider two more measurements. We take into 
account the constraints on VLb from the abundance of light 
elements. We use n B h 2 = 0.022 ± 0.002, where the Hubble 
constant is defined to be Ho = 100/ikm s _1 Mpc _1 . Lastly 
we consider current constraints from the Key Project of the 
Hubble Space Telescope (HST) on the expansion rate today. 
We use H = 72 ± 8km s^Mpc" 1 . 



III. EXPLORING PARAMETER SPACE 

We are interested in seeing if the presence of massive neu- 
trinos is a generic feature of the class of models that we are 
considering. We assume three families of neutrinos with iden- 
tical masses, and we take the mass of each family, m„, as the 
free parameter. An obvious first case to study is a generaliza- 
tion of the TeVeS result from Skordis et al, i.e. a Euclidean 
Universe with a cosmological constant and a constant effec- 
tive Newton's constant. Indeed we find that the posterior for 
m v is positive, centered at m„ ~ 0.84 eV and we can set a 
lower bound on the neutrino mass at the 95% confidence level 
(CL) of m v > 0.31 eV, as shown in table (Q). The value of 



m v proposed in Skordis et al lies comfortably in that range 
0.31<m„<1.48eV. 

Relaxing the assumption that the acceleration is driven by 
a cosmological constant (i.e. freeing up wo and w\), leads 
to a lower bound of m„ > 0.92 eV, slightly stronger than in 
the previous case. The supernovae data strongly constrain the 
parameters describing the nature of dark energy in the range 
< z < 1.8, the era where its contribution is dominant, and 
favour an effective w(z) < —1. This means that the contri- 
bution of the dark energy component diminishes more rapidly 
than A as a function of z. To compensate, the neutrinos are 
required to be relativistic at the surface of last scattering and 
hence considerably more massive, giving rise to the observed 
shift in the distribution to higher masses. The increased free- 
dom gives a broader distribution. When the supernovae data 
set is removed, a larger contribution to the total energy from 
the dark energy component is allowed, weakening the lower 
bound on m„, as shown in Table (TO . However, for this simple 
class of theories, strong statements can now be made: there is 
a definite lower bound on the mass of the neutrino, as can be 
seen from Table 

Up until now, studies of theories of modified gravity have 
been undertaken in the context of Euclidean Universes. Re- 
laxing the assumption of spatial flatness greatly broadens the 
posterior distribution of m„, in particular, extending it to as 
much as 1.8 eV at the 95% CL. These models correspond to 
closed Universes where Qk < 0. In addition the dark en- 
ergy parameters are less strongly peaked (due to the degen- 
eracy with curvature). Models in which the Universe is open 
are however favoured, leading to a weakened lower bound of 
m ll > 0.07 eV but only at the 68% CL. In the absence of the 
supernovae data, Q,k is weakly constrained accordance with 
the increased freedom. This leads to a generally broader m v 
distribution with a similar peak. 

We have been exploring the effect of the extra degrees of 
freedom but we should expect modifications to the left hand 
side of equation Q] We have parameterized this in terms of 
G e f f that depends on the scale factor. In principle, the time 
dependence of G e / / can be more complex, depending on the 
normal matter fields as well as the extra degrees of freedom. 
Furthermore, for any given theory of modified gravity, the 
Bianchi identities as well as the various couplings between 
G e f f and the remaining sector, impose specific constraints on 
its time evolution ll22ll . I.e. we do no have complete freedom 
to vary G eff . 

In what follows, we will be conservative and jettison any 
constraints that come from consistency but we will consider 
two types of relatively general behaviour which encompasse 
what we have found in a wide range of models. One simple 
parametrization is 

G eff = G (l + z) n 

For example for TeVeS one finds that, for a sufficiently small 
n one can adequately mimic the behaviour of G e //. Note that 
this parametrization does lead to a monotonically changing 
G e ff, all the way back to recombination and so must really 
only be considered an approximation- if not, it might lead to 
substantial changes to the peak structure in the CMB at re- 
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combination and we have argued that this is clearly not the 
case. 

A variable G e // can have a substantial effect on allowed 
neutrino masses. In particular, in a Euclidean Universe with 
cosmological constant, it lowers the required mass contribu- 
tion from neutrinos significantly to m v < 0.35 eV at the 
2a level such that the massless case is no longer ruled out. 
This implies an anti-correlation between n and m v . Extend- 
ing the model further to include dark energy again requires 
larger neutrino masses (lower bound of 0.22 eV at 95% CL). 
However in non-flat Universe case, the freedom in parameter 
space means that a wide range of masses are tolerable, includ- 
ing the massless scenario and m„ = 1.32 eV at the 2a level. 
We note that allowing for the possibility of a time-dependent 
G e ff(z) parameterized as above and admitting spatial curva- 
ture will have similar effects on the Hubble equation. The pri- 
mary effect of G e f / is to shift the distribution of m„ to lower 
values, while fix increases the range of neutrino masses that 
can be tolerated. This is explicitly illustrated in Figure (Q]). 
The plot (a) compares the 68% and 95% confidence intervals 
in (m„,0^r) space when G e // is time-independent (dashed 
lines) and dynamical (solid lines) and shows the shift in the 
allowed regions to more negative values of £Ik- Figure (QJ) 
illustrates the impact of spatial curvature on constraints on m„ 
in the presence of a time-dependent G e //. The la and 2a re- 
gions are significantly reduced when curvature is admitted. 

Another possible parameterization is if G e f / switches be- 
tween two values at some point in the past. For example, 
if G e ff is approximately six times larger during the baryon 
dominated era than it is now, the background evolution will 
be essentially equivalent to that of dark matter dominated Uni- 
verse at that time. To mimic this effect we consider 



G 



eff 



G 1 



1 



- 7Z 



(5) 



At low redshift, G e f / starts at Go and increases linearly with 
z. At large redshift (z > I/7), G(z) tends towards a constant 
Go a/7. We limit the change in G(z) from z = to recombi- 
nation by imposing the condition that a/7 < 5 such that its 
does not change by factor of more than 6. We find that with 
this parametrization, which is reminiscent of a number of dif- 
ferent models, that the results are almost identical to that of 
the previous proposal for G e //. Indeed, it is the very late time 
behaviour of G e / / that plays a significant role in changing the 
observables and in that respects the two parametrizations are 
very similar. 



Model 


68% CL 


95% CL 


ACDM 


m v < 0.16 


m v < 0.38 


wCDM 


0.32 < m v < 0.97 


0.064 < m v < 1.27 


iuCDM + n K 


0.02 < m v < 1.27 


m„ < 3.065 


ACDM +G(z) 


m u < 0.14 


m u < 0.35 


toCDM + G(z) 


0.015 < m„ < 0.53 


m v < 0.92 


wCDM +Q K + G(z) 


m u < 1.38 


m v < 2.79 


ACDM +G(a,f,z) 


m v < 0.16 


m v < 0.35 



TABLE II: Results for different cosmological models for a compila- 
tion data set where the supernovae data is excluded. 
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IV. DISCUSSION 

It has been claimed that modified theories of gravity in- 
evitably require the presence of massive neutrinos and that 
these may be sufficiently massive to be measurable with up 
and coming neutrino experiments such as KATRIN ll23ll . This 
claim has been triggered by two pieces of anecdotal evidence. 
Firstly that the simplest TeVeS model needs neutrinos to fit the 
angular power spectrum of the CMB as shown in [10]. And 
secondly, that attempts at reconciling observed and inferred 



FIG. 1: (a) The 68% and 95% confidence regions in (m„, Slpc ) space 
from the compilation data set excluding the supernovae data. The re- 
gion with dashed lines corresponds to a model where n — (ie. 
G(z) = Go), while the region with solid lines corresponds to a Uni- 
verse where n (ie. G(z) = G{z,n)). (b) The 68% and 95% 
confidence regions in (m„ , n) space. The region with dashed lines 
corresponds to a model where Q.k = 0, while the region with solid 
lines corresponds to a Universe where Q.k 7^ 0. The plots illustrate 
the general trend of G e // to shift the distribution of m v to lower 
values, while Q.k increases the range of neutrino masses that can be 
tolerated. 
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masses of clusters requires the presence of a massive neutrino 
halo i24l l25ll . In this letter we have attempted to extend the 
remit of the first piece of evidence. 

We have found that, although for a restricted set of models, 
we can place a lower bound on the mass of the neutrino, for 
more general ranges of parameters, it is possible to satisfy the 
subset of cosmological constraints without having to invoke 
massive neutrinos. This is not to say that specific models with, 
for example, a variable effective Newton's constant might not 
lead to a tight constraint on the neutrino mass. But it is clearly 
not possible to make a definitive statement on the mass of the 
neutrino for general theories of modified gravity. Theories 
must be studied case by case and we have shown how this can 
be done in an economical way. 

It may be possible to come up with constraints on the neu- 
trino masses from a different set of observables, related to the 



second piece of evidence. For example, in the simplest pic- 
ture of a cluster in these theories, neutrinos seem to be in- 
evitable to be able to make up dynamical mass measurements 
and weak lensing observations. This simple picture is incom- 
plete and much of the work that has been done on clusters in 
the context of modified gravity has opted to ignore the extra 
degrees of freedom B26I1 . They can play a significant role and, 
in the same way as for large scale observations, may substan- 
tially weaken cluster constraints on the neutrino mass. A more 
detailed analysis of these systems must be undertaken before 
definitive conclusions can be inferred. 
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